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We study the generation of magnetic fields during preheating within an scenario of hybrid inflation 
at the electroweak (EW) scale. We find that the non-perturbative and strongly out-of-equilibrium 
process of magnetic field production occurs along the lines predicted by Vachaspati many years 
ago. The system starts in the false vacuum at the end of inflation, and very quickly the initial 
quantum fluctuations of the Higgs field get amplified via long wavelength spinodal instabilities. The 
subsequent nucleation of the random Gaussian Higgs field bubbles (lumps) leads to EW symmetry 
breaking, and to the creation of Z-strings, which soon decay, along with longwave magnetic flux 
tubes with nontrivial helicity. The intensity and scales in these helical magnetic fields are consistent 
with their later development into the microgauss fields observed in galaxies and clusters of galaxies. 



The origin of primordial magnetic fields (PMF) is one 
of the remaining puzzles in cosmology [l| . Magnetic fields 
play an important role in the evolution of the primordial 
plasma in the early universe, and in the propagation of 
high energy cosmic rays in our galaxy and beyond. They 
may influence galaxy and large scale structure formation, 
and may generate a stochastic background of gravita- 
tional waves. 

Magnetic fields have been measured on the scales of 
galaxies and clusters of galaxies with a magnitude of or- 
der microgauss [2] . There is even some evidence of their 
presence on scales of superclusters, and associated with 
quasars at redshifts z ~ 2. The main difficulty in under- 
standing their origin is not in their magnitude but in its 
correlation scale. The fiG order of magnitude of galactic 
MF could be explained by an amplification via a dynamo 
mechanism initiated by a tiny seed, B ~ 10 -30 G, in the 
context of gravitational collapse in a flat ACDM model. 
Their magnitude on clusters scales is much more difficult 
to explain via a dynamo mechanism In any case, an 
initial seed will be required. There are both astrophysi- 
cal and cosmological theoretical models trying to account 
for the origin of the primordial seed: Biermann battery 
in intergalactic shocks, stellar magnetic winds (like in our 
Sun), supernova explosions, galactic outflows in the inter- 
galactic medium (IGM), quasar outflows of magnetized 
plasma into the intra-cluster medium (ICM) [1]; as well as 
from early universe phase transitions 0], from magnetic 
helicity at the EW transition, together with the baryon 
asymmetry of the universe (BAU) [4], via hypercharge 
and hypermagnetic field generation before the EW tran- 
sition p], from second order cosmological perturbations 
from inflation [6], from reheating after inflation [7[, or 
produced during inflation [8], etc. 

The ubiquity of microgauss MF on all scales sug- 
gests a cosmic origin. Could they be generated in the 
early universe and then redshifted until today? Mag- 
netic fields could have appeared as the electroweak force 
broke into weak interactions plus electromagnetism, at a 
scale Tew ~ 100 GeV, when the universe had an energy 
density pew ~ 10 8 GeV 4 . The universe was (or became) 
radiation dominated at that time, and the energy density 
carried by magnetic fields could correspond to a fraction 



/ of the total. If we redshift this energy density until to- 
day (T = 2.725 K) we get p B (t ) = (T /T EW ) 4 / Pew ~ 
3 x 10- 53 / GeV 4 = 0.4 /eV/cm 3 = (5 / 1/2 /iG) 2 /(87r), 
where we have used 1 G = 1.95 x 10 -20 GeV 2 . This would 
be enough to explain the cluster and supercluster values, 
and would perhaps require a mild dynamo mechanism to 
amplify it to galactic values (if the fraction / < 1). 

However, a priori, it is not so clear how one can ob- 
tain the large correlation length of magnetic fields ob- 
served at galactic and cluster scales. Any physical mech- 
anism that creates magnetic fields must necessarily be 
causal, but at high temperatures in the early universe 
there is also a natural coherence scale given by the phys- 
ical horizon. At the EW scale the particle horizon is 
10 -10 lightseconds (~ 3 cm), which today corresponds 
to a comoving scale of ~ 1 AU, clearly insufficient when 
compared even with the irregular (turbulent) component 
of the galactic magnetic field (L ~ 100 parsecs), not to 
mention the regular (uniform) component, which has cor- 
relations L ~ 10 kpc. It thus seems apparently impos- 
sible to explain with our mechanism the coherent mag- 
netic fields observed on galaxy clusters and superclus- 
ter scales (of order 10 Mpc) with intensities of order /aG 
to nG. Nevertheless, if we assume that the plasma after 
the EW transition is sufficiently turbulent to maintain 
magnetic fields of the largest possible coherence scales 
via inverse cascade, then we could reach cosmological 
scales today. The largest coherence scale is the physi- 
cal horizon. If a strong inverse cascade is active, then 
the coherence length of the magnetic fields could grow 
as fast as the horizon, i.e. like the scale factor squared 
during the radiation dominated era. This optimal sit- 
uation is only attainable in the presence of a plasma, 
and thus is bound to stop acting at photon decoupling, 
when the universe becomes neutral. Since then, the cor- 
relation length can only grow with the expansion of the 
universe, as the scale factor. Taking these facts into ac- 
count, and using the adiabatic expansion relation, we 
can compute the coherence scale of the field in terms of 



that at the electroweak scale, £o = £ew (f^w) 

6 x 10 25 cm = 20 Mpc, where we have 
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made the approximation that matter-radiation equality 
and photon decoupling occurred more or less simultane- 
ously (a careful computation gives only a minor correc- 
tion). The surprising thing is that this simple calculation 
gives precisely the order of magnitude for the largest cor- 
relation length of cosmic MF ever observed (i.e. cluster 
scales). If the jiG magnitude of the primordial MF seed, 
arising from simple scaling down the energy density since 
the EW transition, seemed a curious coincidence, the fact 
that an inverse cascade could also be responsible for the 
observed correlation length becomes a suspicious coinci- 
dence. These observations have triggered a large number 
of investigations devoted to the analysis of magnetogen- 
esis during the electroweak phase transition 0,13, [l| . 

In this Letter we address the question of whether the 
required conditions -a significant fraction of energy den- 
sity in magnetic fields with a large helical component 
and sustained periods of inverse cascade - can be ob- 
tained within a simple and concrete scenario of low-scale 
hybrid inflation [9[ . We analyze whether large scale mag- 
netic fields are generated during the preheating and re- 
heating periods which lead to EW symmetry breaking. 
Previous attempts to analyze magnetogenesis during re- 
heating as a consequence of spinodal instability can be 
found in Ref. [7] . Our approach differs significantly and is 
based on the use of the classical approximation: integra- 
tion of the non-linear equations of motion with stochastic 
initial conditions determined by matching the quantum 
expectation values at the end of inflation. The numerical 
implementation is based upon discretization of our sys- 
tem onto a space-time lattice. Our results indicate that a 
mechanism resembling that introduced by Vachaspati [4] 
could have taken place in the early universe. The mech- 
anism relies on the fact that an inhomogeneous Higgs 
background may act as a source for magnetic fields. The 
effect would come from the existence of non-trivial gra- 
dients of the Higgs field phase which would generate W- 
currents and magnetic fields during a cold EW transi- 
tion Q. These gradients are usually associated with the 
appearance of unstable Z-strings, together with the he- 
lical magnetic field necessary for growth into large scales 
via MHD turbulence [4j]. As we will see, our results show 
that during the initial stage of electroweak symmetry 
breaking (EWSB), our setup is indeed characterized by 
the appearance of Z-strings and magnetic strings with 
non-trivial helicity. 

We consider a specific model of Hybrid inflation with 
the field content of the gauge plus scalar sector of the 
Standard Model coupled to a singlet scalar field, the in- 
flaton The gauge sector includes both the SU(2) and 
the hypercharge U(l) fields with corresponding gauge 
couplings gw and gy- The scalar sector contains the 
4-real component Higgs field M (x), and a singlet infla- 
ton field x( x ) coupled to the Higgs via the term g 2 ^ x 2 • 
Our initial condition assumes that at the end of infla- 
tion both the Higgs and the gauge fields are in the de 
Sitter vacuum. In fact, since for EW-scale inflation the 
rate of expansion is negligible compared to any other 




FIG. 1: The location of the Higgs lumps (light/red) and the 
magnetic field flux tubes (dark/blue) at mt =15. 

scale, H « 1(T 5 eV <C v = 246 GeV, the de Sitter 
vacuum state is indistinguishable from the Minkowski 
vacuum for the range of momenta we are considering. 
The inflaton field is given by the homogeneous mode 
alone, which is decreasing with time. This time de- 
pendence translates, through the Higgs-inflaton coupling, 
into a time dependence of the Higgs mass-squared, which 
changes from positive to negative inducing an exponen- 
tial growth of the Higgs low momentum modes in a pro- 
cess known as "tachyonic preheating" [10]. The remain- 
ing non-tachyonic fields are populated through their in- 
teraction to the Higgs. Fermionic fields can be safely 
ignored from the analysis since the perturbative decay of 
the Higgs into quarks and leptons occurs within a time 
scale mt ~ 10 3 — 10 4 , many times bigger than the scales 
we are considering in our simulations. 

There have been several numerical studies of this 
model [TO, fllL in an attempt to explain the BAU in 
this contextjO, EH, 0, as well as the production of 
a gravitational wave background [151 ] . In this paper, 
the Hypercharge gauge field is included for the first 
time (preliminary results were discussed in Ref. [13] ) . This 
will allow us to study the generation of U(l) electromag- 
netic fields during preheating. In order to analyze the 
results we first have to provide a definition of the U(l) em 
content of the SU(2)xU(l) fields in the Lagrangian. We 
will be using a gauge invariant definition of Z vector po- 
tential on the lattice as 



4w 



Tr 



2|0(n)||0(n + /i)| 



(1) 



where 4>^{n) are the lattice Higgs field components, with 
modulus |0(77-) I . In matrix notation, <£(n) = (j)oI + i(j) a T a . 
We have introduced the adjoint unit vector n = n a r a , 
with components n a (n) = ip^ (n)r a (p(n) /\(p(n)\ 2 , with 
tp(n) = <I>(n)(l,0) T the Higgs doublet, and is the 
SU(2)xU(l) lattice covariant derivative. Note that our 
definition of the Z vector potential corresponds to the 
standard one in the unitary gauge. We can then com- 
pute the field strength of the U(l) em field in terms of the 
Z and Hypercharge field strengths: 



F;» = S m 2 (9 w )F^(n) - f£ (n) . 



(2) 
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FIG. 2: 2D contour plots of the Higgs norm (top left), mag- 
netic field flux tubes (top right), Z-string helicity (bottom 
left) and magnetic field helicity (bottom right) at mt = 15. 
Note the precise spatial correlation. 

This procedure preserves the U(l) em gauge-invariance on 
the lattice. 

As mentioned previously, there are two essential ingre- 
dients in Vachaspati's proposal for magnetogenesis: the 
existence of inhomogeneities in the Higgs field and the 
generation of a non-trivial magnetic field helicity. Both 
of them are present in the first stages of evolution right 
after inflation ends. This might seem surprising since 
initially both hypercharge and SU(2) magnetic fields are 
zero. The way it comes about is by noticing that local 
variations of the orientation of the Higgs field induce via 
Eq. ([2]) a non- vanishing contribution to the electromag- 
netic field. The initial conditions for the classical evo- 
lution are given by a Higgs field as a Gaussian random 
field. This initial distribution carries essential informa- 
tion about how the process of EW symmetry breaking 
takes place. It is, for instance, the seed for the growth of 
bubble shells which arise from lumps in the initial distri- 
bution of the Higgs field norm. We have observed that 
it also determines the way in which magnetic fields are 
seeded. At later times, after EWSB has occured, the 
correlation between the Higgs field norm and magnetic 
fields is manifest. In Fig. [1] we show a snapshot of both 
quantities at a time (mt = 15) after EWSB has occurred 
and bubble shells (in red) fill almost all the volume of the 
box. Magnetic fields (in blue) are localized in string-like 
structures concentrated in the region between bubbles. It 
is there where the Higgs field remains closer to the false 
vacuum for a longer period of time giving rise to larger 
gradients in the Higgs field phases. This linkage between 
magnetic strings and Higgs field minima is even more ev- 
ident in the two dimensional contour plots presented in 
Fig. [2j As mentioned previously part of this structure 
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FIG. 3: Time evolution of the Higgs v.e.v. (black) and the 
square of the helicity in Z (green) and U(l) (red) magnetic 
fields normalized to the physical volume V. The Z-helicity has 
been rescaled by tan 2 0w> For rah = 2ra w , the corresponding 
slopes are: Photon field: 0.68(7) and Z field: -0.27(4). 

was already present in our initial condition (mt = 5), in 
the absence of SU(2) and hypercharge magnetic fields, 
but it has become far more complex at this stage. 

It is interesting to point out that previous analysis of 
the cold EW transition, neglecting the effect of hyper- 
charge, have revealed that regions between bubbles are 
also sources of Chern-Simons number through the ap- 
pearance of sphaleron-like configurations attached to the 
location of zeroes in the Higgs field norm [H, 03 • Since 
for non-zero Weinberg angle sphalerons are like magnetic 
dipoles, it is tempting to correlate the observed U(l) 
magnetic strings to the alignment of sphaleron dipoles. A 
detailed investigation of this issue is beyond the scope of 
this paper. In any case, the connection between Chern- 
Simons number and magnetic structures is more explicit 
when analyzing the generation of magnetic helicity, which 
can be regarded as the Abelian counterpart of the Chern- 
Simons number. During these first stages of evolution 
the magnetic and Z-boson helicities are indeed tracking 
the magnetic strings (a clear example is shown in Fig. [21 
comparing top and bottom boxes). The proportionality 
of the Z-boson and photon helicities in the initial condi- 
tion is a direct consequence of our definition of the fields. 
Nonetheless this correlation is still preserved once gauge 
fields and non-linearities have started to play a role. 

The previous figures correspond to a snapshot at mt = 
15, when the string structure comes out to be stronger. 
It is important to ascertain whether the later stages of 
evolution preserve relic magnetic fields with non trivial 
magnetic helicity. This quantity is defined as: 

H = [ d?xA.B = i J2 l£ ' ^(V x B*(k)) , (3) 
J k \ k \ 

where V is the volume of space. Being a pseudoscalar and 
since we have not included CP violation, the expectation 
value of the helicity is zero. Nevertheless, to have an 
idea of the efficiency of this mechanism in generating non- 
trivial magnetic helicity we have computed its dispersion. 
Its time evolution is presented in Figj3] where we also 
present the time dependence of the helicity dispersion in 
the Z-boson magnetic field. The latter has been rescaled 
by tan 2 0w to make it agree with the electromagnetic 
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helicity for the initial configuration at mt = 5. While 
initially both helicities coincide, time evolution damps 
the Z-helicity while it enhances the Abelian magnetic 
one. This agrees with Vachaspati's picture, that pre- 
dicted an initial correlation among both quantities, and 
subsequent decay of Z-strings, leaving remnant magnetic 
fields with non-trivial helicity. From our analysis we de- 
termine that the magnetic helicity dispersion is actually 
growing in time with a power-law behavior, ~ ^0.68(7) 
A detailed study of the late time behavior of the system 
will be reported in Ref. [l6[. We will discuss there the 
nature of the plasma during the first stages of reheating 
which shows evidence of electric charge separation. We 
will also present some evidence of the presence of inverse 
cascade in the magnetic power spectrum and the helicity 
dispersion. 

In summary, in this paper we studied the preheating 
and early reheating periods of a model of hybrid infla- 
tion at the EW scale within the classical approximation. 



This is the first study of this sort including SU(2) xU(l). 
We have focused upon the generation of magnetic fields 
in this context and provided evidence that the mecha- 
nism proposed in 1991 by Vachaspati [4| is at work. Spa- 
tial gradients of the Higgs field, seeded by those already 
present in the random Gaussian initial condition, act as 
sources of the magnetic field. The latter has a non-trivial 
helicity dispersion, correlated initially with the presence 
of Z-strings, that is growing with time. 

Acknowledgements - It is a pleasure to thank A. 
Achucarro, M. Hindmarsh, M. Salle, M. Shaposhnikov, 
J. Smit, J. Urrestilla and T. Vachaspati for generous 
discussions. We acknowledge financial support from 
the Madrid Regional Government (CAM) under the 
program HEPHACOS P-ESP-00346, and the Spanish 
M.E.C. under contracts FPA2006-05807, FPA2006-05485 
and FPA2006-05423. We also acknowledge the use of the 
MareNostrum supercomputer and the IFT computation 
cluster. 



[1] E. N. Parker, "Cosmological Magnetic Fields," 
Clarendon Press, Oxford (1979); A. A. Ruzmaikin, 
A. M. Shukurov and D. D. Sokoloff, "Magnetic Fields 
of Galaxies," Kluwer Academic Publisher, Dordrecht 
(1988); P. P. Kronberg, Rept. Prog. Phys. 57, 325 
(1994); K. Enqvist, Int. J. Mod. Phys. D 7, 331 (1998); 

D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 
(2001); L. M. Widrow, Rev. Mod. Phys. 74, 775 (2003); 
M. Giovannini, Int. J. Mod. Phys. D 13, 391 (2004); 
M. Shaposhnikov, AIP Conf. Proc. 784, 423 (2005); 
R. Durrer, New Astron. Rev. 51, 275 (2007). 

[2] J. Vallee, Fund. Cosm. Phys. 19, 1 (1997). 

[3] C. J. Hogan, Phys. Rev. Lett. 51, 1488 (1983); 
J. M. Quashnock, A. Loeb and D. N. Spergel, Astro- 
phys. J. 344, L49 (1989). B. L. Cheng and A. V. Olinto, 
Phys. Rev. D 50, 2421 (1994). G. Baym, D. Bodeker and 
L. D. McLerran, Phys. Rev. D 53, 662 (1996) A. Bran- 
denburg, K. Enqvist and P. Olesen, Phys. Rev. D 54, 
1291 (1996) G. Sigl, A. V. Olinto and K. Jedamzik, Phys. 
Rev. D 55, 4582 (1997); J. Ahonen and K. Enqvist, Phys. 
Rev. D 57, 664 (1998); M. Hindmarsh and A. Everett, 
Phys. Rev. D 58, 103505 (1998); D. T. Son, Phys. Rev. D 
59, 063008 (1999); M. Christensson, M. Hindmarsh and 
A. Brandenburg, Astron. Nachr. 326, 393 (2005). 

[4] T. Vachaspati, Phys. Lett. B 265, 258 (1991); T. W. B. 
Kibble and A. Vilenkin, Phys. Rev. D 52, 679 (1995); 

E. J. Copeland and P. M. Saffin, Phys. Rev. D 54, 6088 
(1996); D. Grasso and A. Riotto, Phys. Lett. B 418, 258 
(1998); E. J. Copeland, P. M. Saffin and O. Tornkvist, 
Phys. Rev. D 61, 105005 (2000); T. Vachaspati, Phys. 
Rev. Lett. 87, 251302 (2001) 

[5] M. Joyce and M. E. Shaposhnikov, Phys. Rev. Lett. 
79, 1193 (1997); M. Giovannini and M. E. Shaposh- 
nikov, Phys. Rev. D 57, 2186 (1998); M. Giovannini and 
M. E. Shaposhnikov, Phys. Rev. D 62, 103512 (2000). 

[6] K. Takahashi, K. Ichiki, H. Ohno and H. Hanayama, 
Phys. Rev. Lett. 95, 121301 (2005); K. Ichiki, K. Taka- 
hashi, H. Ohno, H. Hanayama and N. Sugiyama, Science 
311, 827 (2006) 



[7] E. A. Calzetta and A. Kandus, Phys. Rev. D 65, 063004 

(2002) ; D. Boyanovsky, M. Simionato and H. J. de Vega, 
Phys. Rev. D 67, 023502 (2003); D. Boyanovsky, H. J. de 
Vega and M. Simionato, Phys. Rev. D 67, 123505 (2003); 
D. Boyanovsky and H. J. de Vega, AIP Conf. Proc. 784, 
434 (2005). 

[8] M. S. Turner and L. M. Widrow, Phys. Rev. D 37, 2743 
(1988); M. Giovannini and M. E. Shaposhnikov, Phys. 
Rev. D 62, 103512 (2000). 
[9] J. Garcia- Bellido, D. Y. Grigoriev, A. Kusenko and 
M. E. Shaposhnikov, Phys. Rev. D 60, 123504 (1999). 

[10] G. Felder, J. Garcfa-Bellido, P. Greene, L. Kofman, A. 
D. Linde and I. Tkachev, Phys. Rev. Lett. 87, 011601 
(2001); G. Felder, L. Kofman and A. D. Linde, Phys. 
Rev. D 64, 123517 (2001). 

[11] J. Garcfa-Bellido, M. Garcia Perez and A. Gonzalez- 
Arroyo, Phys. Rev. D 67, 103501 (2003). 

[12] J. Garcfa-Bellido, M. Garcia Perez and A. Gonzalez- 
Arroyo, Phys. Rev. D 69, 023504 (2004). 

[13] J.Smit, A.Tranberg, JHEP 0212, 020 (2002); JHEP 
0311, 016 (2003); JHEP 0608, 012 (2006). J. I. 
Skullerud, J. Smit and A. Tranberg, JHEP 0308, 045 

(2003) ; M. van der Meulen, D. Sexty, J. Smit and 
A. Tranberg, JHEP 0602, 029 (2006); A. Tranberg, 
J. Smit and M. Hindmarsh, JHEP 0701, 034 (2007); 
J. Smit, in proceedings of Lattice 2005 conference, PoS 
(LAT2005) 022, and references therein. 

[14] A. Dfaz-Gil, J. Garcfa-Bellido, M. Garcia Perez, 
A. Gonzalez- Arroyo, in proceedings of Lattice 2005 
conference, Po S (LAT2005) 242 (2006) [arXiv: 
hep-lat/0509094]; ibid, in proceedings of Lattice 2007 
conference, PoS (LAT2007) 052 f2007). larXiv:0710.0580l 
[hep-lat] . 

[15] J. Garcia-Bellido and D. G. Figueroa, Phys. Rev. Lett. 

98 (2007) 061302; J. Garcia -Bellido, D. G. Figueroa and 

A. Sastre, larXiv: 0707. 08391 [hep-ph] . 
[16] A. Dfaz-Gil, J. Garcfa-Bellido, M. Garcia Perez, A. 

Gonzalez- Arroyo, in preparation. 



